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Abstract 

| The topological defect scenario of origin of the observed highest energy cosmic rays is 

reviewed. Under a variety of circumstances, topological defects formed in the early Universe 
can be sources of very massive particles in the Universe today. The decay products of these 
massive particles may be responsible for the observed highest energy cosmic ray particles 
above 10 20 eV. Some massive particle production processes involving cosmic strings and 
magnetic monopoles are discussed. We also discuss the implications of results of certain 
recent numerical simulations of evolution of cosmic strings. These results (which remain to 
be confirmed by independent simulations) seem to show that massive particle production 
may be a generic feature of cosmic strings, which would make cosmic strings an inevitable 
source of extremely high energy cosmic rays with potentially detectable flux. At the same 
time, cosmic strings are severely constrained by the observed cosmic ray flux above 10 20 eV, 
if massive particle radiation is the dominant energy loss mechanism for cosmic strings. 

Introduction 

Cosmic Topological Defects (TDs)fl], § such as magnetic monopoles, cosmic strings, domain 
walls, and various hybrid TD systems consisting of these basic kinds of TDs, are predicted 
to form in the early Universe as a result of symmetry-breaking phase transitions envisaged 
in unified theories of elementary particle interactions. Under a variety of circumstances 
these TDs can be sources of extremely massive unstable particles in the universe today^, 
f|, |^, H [7|, |^, [It], |T1|] . The masses m x of these so-called "X" particles (the quanta of 



1 Invited talk given at the Workshop on "Observing the Highest Energy Particles (> 10 20 eV) from Space" , 
College Park, Maryland, USA, November 13 - 15, 1997, to be published in the Proceedings (AIP). 

2 NAS/NRC Resident Senior Research Associate at NASA/GSFC on sabbatical leave from Indian Institute 
of Astrophysics, Bangalore-560 034. India. 



1 



the massive gauge- and higgs fields of the underlying spontaneously broken gauge theory) 
would typically be of order the symmetry-breaking energy scale at which the relevant TDs 
were formed, which, in Grand Unified Theories (GUTs), can be as large as ~ 10 16 GeV. 
The decay of these X particles can give rise to extremely energetic photons, neutrinos and 
nucleons with energy up to ~ nix- If the X particle production rate from TDs is large 
enough, these extremely energetic particles may be detectable by ground-based as well as 
space-based large-area detectors being planned for detecting ultrahigh energy (UHE) (i.e., 
energy > 10 9 GeV) cosmic rays. These cosmic ray detectors may thus provide us with a tool 
for studying the signature of TDs and thus of GUT scale physics. 

There is currently much interest in the possibility that the Extremely High Energy (EHE) 
Cosmic Ray (EHECR) events — those with energies above 10 11 GeV reported recently ||12|| 
- may be due to decays of massive X particles originating from TDs[f[3|, [II], [14], [L5|, [16], [T7|, 



1 1 . This possibility is of interest in view of the fact that the energies associated with the 



EHECR events are hard to obtainflTSl O, P1I within the standard diffusive shock acceleration 



mechanism [jJIJ that involves first-order Fermi acceleration of charged particles at relativistic 
shocks associated with known powerful astrophysical objects; see, however, Ref. [ETJ. In 



addition, there is the problem of absence of any obviously identifiable sources for these 



EHECR events p2|, [13]. These problems are avoided in the TD scenario in a natural way. 
Firstly, no acceleration mechanism is needed: The decay products of the X particles have 
energies up to ~ mx which can be as large as, say, 10 16 GeV. Secondly, the absence of 
obviously identifiable sources is not a problem because TDs need not necessarily be associated 
with any visible or otherwise active astrophysical objects such as AGNs or radio galaxies. 

The basic ideas of the TD scenario of origin of EHECR have been reviewed in a number 
of discussions in the past; see, e.g., Refs. |F| ^3], [24], p5| . Detailed calculations of the 



predicted spectra of nucleons, photons, and neutrinos in the TD scenario have been done 
in the past several years|^, |26], [10], [15|, [H], [T7|, [16], |27|. Constraints on the TD scenario 



imposed by experimental data on EHECR and diffuse gamma ray background have also 
been discussed [[28], |lj, |T7|, [IjJ. 

In what follows, I first discuss briefly some of the basic aspects of topological defects and 
their formation in the early Universe. The basic steps in the calculation of the 'observable' 
particle spectra resulting from the decay of the X particles is discussed. A simple benchmark 
calculation of the X particle production rate required to explain the observed EHECR particle 
flux is performed. I then discuss three specific X particle production mechanisms involving 
(1) "ordinary" cosmic strings, (2) monopolonia — metastable monopole-antimonopole bound 
states, and (3) cosmic "necklace" — a system of monopoles on strings. I also discuss the 
implications of the results of certain recent numerical simulations of evolution of cosmic 
strings in the Universe. These results, if confirmed by independent simulations, would imply 
that massive X-particle production and hence cosmic ray production might be a generic 
feature of cosmic strings, which would make cosmic strings an inevitable source of EHECR 
with potentially detectable flux. Indeed, in this case, we shall see that the measured EHECR 
flux already puts severe constraint on the energy scale of symmetry-breaking associated with 
any cosmic string forming phase transition in the early Universe. 
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I use natural units with K — c — ks — 1 throughout, unless otherwise stated. 



Topological Defects and X Particle Production: A Brief 
Review 

Topological Defects are sometimes characterized as "exotic". In actual fact, TDs are rou- 
tinely seen, measured, and studied in condensed matter systems in laboratories. TDs form 
during phase transitions associated with the phenomenon of spontaneous symmetry breaking 
(SSB), which is a central concept in condensed matter physics as well as in the Standard 
Model of particle physics. Well-known examples of TDs in condensed matter systems are 
vortex lines in superffuid helium, magnetic flux tubes in type-II superconductors, disclina- 
tion lines and 'hedgehogs' in nematic liquid crystals, and so on. Perhaps what is perceived 
as exotic is the existence of TDs in the cosmo logical context. However, it has been real- 
ized for quite some time now, particularly since the early seventies, that our Universe in its 
early stages must have behaved very much like condensed matter systems. Indeed, ideas of 
unified gauge theories of elementary particle interactions taken together with the hot big- 
bang model of the early Universe necessarily imply that our Universe in its early history has 
passed through a sequence of symmetry-breaking phase transitions as it expanded and cooled 
through certain critical temperatures. Depending on the symmetry breaking pattern, one 
or more kinds of the three basic kinds of topological defects — magnetic monopoles, cosmic 
strings and domain walls — could be formed during some of these phase transitions P], H. In 
fact, formation of magnetic monopoles, is inevitable in practically all Grand Unified Theories 
(GUTs) that provide a unified description of the electroweak and strong interactions^. The 
monopoles are analogous to the 'hedgehogs' in nematic liquid crystals and appear whenever 
the unbroken symmetry group possesses a local U(l) symmetry. The 'global' cosmic strings, 
which arise in breaking of global U(l) symmetry, are similar to vortex filaments in super- 
fluid helium, and the 'local' or 'gauge' cosmic strings arising from breaking of a local U(l) 
symmetry are similar to the magnetic flux tubes in type-II superconductors. Cosmic domain 
walls appear whenever a discrete symmetry is spontaneously broken. 

Interestingly, it has recently become possible to simulate the analogue of cosmic string 



formation in the early Universe by means of laboratory experiments [30] on vortex-filament 
formation in the superfluid transition of 3 He, which occurs at a temperature of a few mil- 
likelvin. The results of these experiments have provided striking confirmation of the basic 
Kibble-Zurek picture |2|, |3l| of topological defect formation in general, which was initially 



developed within the context of defect formation in the early Universe]]]]. 

3 The inevitability of monopole formation leads to the well-known "monopole problem" of cosmology, 
which historically was one of the "problems" that motivated the idea of inflationary cosmology; for a review, 
seeRef.|29tl 
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Symmetry Breaking Phase Transitions and Formation of Topolog- 
ical Defects 

During a symmetry breaking phase transition, the system under consideration undergoes a 
transition from a state of higher symmetry to one of a lower (reduced) symmetry at a critical 
temperature during the cooling of the system. In spontaneous symmetry breaking (SSB), 
the system below the critical temperature possesses multiple degenerate ground states rather 
than a unique ground state. These degenerate ground states differ from each other by the 
'phase angle' or some internal degrees of freedom of the "order parameter" field (OPF) whose 
absolute magnitude (which is same for all the ground states) is a measure of the order (or 
lack of symmetry) in the system. (In particle physics, the OPF is the "higgs" field.) The 
symmetry under consideration is invariance of the energy (or more precisely the Lagrangian 
for the OPF or the higgs field) under transformations that change the 'phase' of the OPF. 
The Lagrangian is always invariant under these transformations because, by construction, it 
depends only on the absolute magnitude of the OPF and not on its 'phase'. However, the 
ground states transform among each other under the action of these transformations, and so 
any chosen ground state is clearly not invariant under the transformations — the symmetry 
is spontaneously broken. The existence of multiple degenerate ground states is the defining 
characteristic of the phenomenon of spontaneous symmetry breaking. By convention, the 
absolute value of the OPF is taken to be zero in the high-temperature unbroken-symmetry 
phase and unity in the low-temperature broken symmetry phase. 

Because of the availability of these multiple degenerate ground states in the low tem- 
perature phase, different parts of the system may choose to settle down to different ground 
states when making the transition to the low temperature phase, especially so if the transi- 
tion happens in an out-of-equilibrium situation. Indeed, one expects that the choice of the 
phase of the OPF in regions separated by more than the correlation length of the thermal 
fluctuations of the OPF will be uncorrelated. The correlation length will always be finite in 
a finite physical system. In the context of the Universe as a whole, there is also an upper 
limit to the correlation length at any time t, namely, the causal horizon length ~ ct. Thus 
the choice of the phase of the OPF will be random and in general different in different parts 
of the Universe separated by more than the causal horizon length at the time of phase tran- 
sition. This often leads to 'obstruction' in the way of uniform completion of the transition 
throughout the bulk of the system. Indeed, it is often the case that the random choice of 
different ground states in different regions leads to some regions being forced to remain in 
the unbroken symmetry phase. These regions are the 'topological defects'. 

The topological nature of the defects becomes clear when one considers the configuration 
of the OPF in the low temperature phase at the end of the phase transition. The choice 
of the 'phase' of the OPF corresponding to different ground states in different parts of the 
system may be such that, in order to avoid energetically unfavored discontinuity in the spatial 
variation of the OPF, the magnitude of the OPF is forced to be zero on some geometrical 
points, lines, or surfaces, which define the 'center' of the defects. Actually, a defect has a 
finite size dictated by the need to minimize the overall energy of the system; the absolute 
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value of the OPF increases gradually from zero at the center to its broken-symmetry-phase 
value at a finite distance from the center. 

The topological stability of the defects is due to non-trivial topological 'winding' of the 
OPF configuration around the defect center. For example, in the case of the linear defect 
(the vortex filament in the superfluid helium, for example), the OPF (the wave function of 
the condensed helium atoms, for example) is a complex number whose phase turns by an 
integral multiple of 2tt as one makes a complete circuit along a closed curve around any 
point on the defect line (the closed curve being on the normal plane cutting the defect line 
at the given point). Once formed, such a configuration cannot 'unwind' by itself and is thus 
topologically stable. 

In the context of spontaneously broken gauge theories, explicit analytical and/or numer- 
ical finite-energy, extended, topologically stable solutions for the higgs- and gauge field con- 
figurations representing cosmic strings, monopoles and domain walls are known; see Ref.@ 
for review. In the broken-symmetry phase, a higgs field responsible for spontaneous symme- 
try breaking of a local gauge theory is massive, as are the gauge bosons of the theory, the 
mass scale being set by the absolute magnitude of the higgs field (more precisely, its vacuum 
expectation value) in the broken symmetry phase. The size of the 'core' of a defect is of 
order m^ 1 , where X represents the higgs or the gauge field. Within the core, the symmetry 
remains unbroken, and the energy density (associated with the higgs and the gauge fields) 
is higher within the core than outside. The topological stability of the defect ensures the 
'trapping' of the excess energy within the core of the defect, which is what makes a defect 
massive. It can be shown that the mass scale of a defect is fixed by the energy (or tempera- 
ture) scale of the symmetry breaking phase transition at which the defect is formed. Thus, if 
we denote by T c the critical temperature of the defect-forming phase transition in the early 
Universe, then the mass of a monopole formed at that phase transition is roughly of order 
T c , the mass per unit length of a cosmic string is of order T c 2 , and the mass per unit area of 
a domain wall is of order T c 3 . 

The X particles are the quanta of excitations of the higgs and gauge fields. In the broken 
symmetry phase, these quanta are massive, their mass is also roughly of order T c . These 
massive quanta typically have very short life times, and so they all decayed away quickly 
soon after the phase transition in the early Universe, and none of those X particles survive in 
the present Universe. This is to say that outside of a defect in the low-energy Universe today, 
the higgs and the gauge fields are in their ground states and no excitations of the X particle 
quanta are present. However, inside the core of a defect, the symmetry is unbroken — the 
X particles are massless inside. Topological stability of the defect prevents this "piece of 
the early Universe" trapped inside the defect from decaying. If, however, there is a process 
which removes this topological protection, then the energy trapped inside the defect will 
dissociate into quanta of the X particles which, being now in the broken-symmetry phase, 
will be massive and short lived, decaying into elementary particles such as quarks and leptons 
which, in turn, would eventually materialize into energetic nucleons, photons, neutrinos, etc., 
that might contribute to the observed EHECR flux. 

Production of X particles from TDs may happen in a variety of ways directly or indi- 
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rectly related to local removal of topological stability of (parts of) TDs. Examples include 
'cusp' evaporation from cosmic strings @, collapse of macroscopic cosmic string loops[|J, 
shrinking of cosmic string loops to radii of order m x l @, annihilation of a monopole with 
an antimonopole[|], [K|, and so on. In the case of current-carrying superconducting cos- 
mic strings|32|, IS], the charge carriers (which could be quanta of a superheavy fermion field 
trapped in 'zero mode' inside the string, or a charged scalar field living inside the string due 
to energetic reasons) are expelled from the string when the current on the string reaches a 
critical value — in the vacuum away from the string, these charge carriers act as the massive 



X particles. In the case of 'ordinary' cosmic strings, recent field theory simulations |33| of 
evolution of cosmic strings in the early Universe show that a cosmic string network loses 
energy directly into oscillations of the underlying gauge and higgs fields 'constituting' the 
string, which in quantum theory, correspond to quanta of massive X particles — a result 
which has important implications for both cosmic strings as well as for EHECR; this result, 
however, remains to be confirmed by independent simulations. 

X-Particle Production Rate 

The number density of X particles produced by TDs per unit time, dnx/dt, can be generally 
written as[||] 

dt w m x \t J y ' 

where t denotes the present epoch, and Q is the rate of energy density injected in the 
form of X particles in the present epoch. The quantity Qo and the parameter p depend on 
the specific TD process under consideration. In writing Eq. (HD it is assumed that the only 
time scale in the problem is the hubble time t and that any other time scale involved can be 
expressed in terms of the hubble time. Similarly, we assume that any energy scale involved 
in the problem is expressible in terms of the energy scale 77 of the symmetry breaking at 
which the TDs under consideration are formed. (Note that nix is fixed by rj.) These 
assumptions are sometimes expressed by saying that the TDs under consideration evolve 
in a scale independent way. Indeed, it turns out that Eq. ([I]) is a phenomenologically 
good parametrization for the specific TD processes studied so far. For example, p = 1 
for a process of X particle production from cosmic string loops [Q, for a process involving 
collapsing monopole-antimonopole bound states ||, |10| , for a process of particle production 
from monopole-string systems called "necklaces" [11]], and so on, while p = for a process 
involving superconducting cosmic string loops ||. 

Note that X particle production from TDs may occur continually at all epochs after 
the formation of the relevant TDs in the early Universe. However, only those X particles 
produced in the relatively recent epochs and at relatively close-by, non-cosmological distances 
(< 100 Mpc) are relevant for the question of EHECR. This is because, nucleons of energies 
above 10 11 GeV produced by the decay of X particles occurring at distances much larger than 
~ 50 Mpc suffer drastic energy loss during their propagation, due to photopion production 
on the the cosmic microwave background radiation (CMBR) fields (the so-called "GZK 
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effect" |[34|), and hence do not survive as EHECR particles today. Distance of sources of 
photons of energies above ~ 10 11 GeV are also similarly restricted due to absorption through 
e + e~ production on the radio background photons (see, e.g., [PS|). The neutrinos, however, 
can survive from much earlier cosmological epochs, and may, if detected, prove to be the 
ultimate discriminant between a TD scenario and a more conventional scenario of origin of 
EHECR. 

From X Particles to 'Observable' Particles 

The X particles released from TDs would decay typically into quarks and leptons. The 
life-time r is typically ~ (arrix)~ 1 (where a ~ few x 10~ 2 ), which for nix ~ 10 16 GeV is 
~ 10 -39 sec or so. The decay is, therefore, essentially instantaneous at late cosmological 
epochs of interest to us. The quarks would hadronize (typically on a strong interaction time 
scale ~ 10 _23 sec, i.e., again practically instantaneously) by producing jets of hadrons, most 
of which would eventually be light mesons (pions) with a small admixture (typically < 10%) 
of baryons and antibaryons (nucleons and antinucleons). The neutral pions decay to two 
photons, while the decay of charged pions gives rise to neutrinos. Some leptons (charged 
as well as neutral) could also be produced directly from the X particle decay. But by far 
the largest number of nucleons, photons, and neutrinos, would be produced through the 
hadronic channel. The spectra of produced particles are, therefore, essentially determined 
by the process of fragmentation of quarks/gluons into hadrons as described by QCD. 

Quark — ► Hadron Fragmentation Spectrum 

The exact process by which a single high energy quark gives rise to a jet of hadrons is 
not known; it involves some kind of non-perturbative physics that is not well understood. 
However, different semi-phenomenological approaches have been developed which describe 
the hadronic "fragmentation spectrum" of quarks/gluons that are in good agreement with 
the currently available experimental data on inclusive hadron spectra in quark/gluon jets in 
a variety of high energy processes. 

In these approaches, the process of production of a jet containing a large number of 
hadrons, by a single high energy quark (or gluon), is 'factorized' into two stages. The first 
stage involves 'hard' processes involving large momentum transfers, whereby the initial high 
energy quark emits 'bremsstrahlung' gluons which themselves create more quarks and gluons 
through various QCD processes. These hard processes are well described by perturbative 
QCD. Thus a single high energy quark gives rise to a 'parton cascade' — a shower of quarks 
and gluons — which, due to the high energy nature of the process, is confined in a narrow 
cone or jet whose axis lies along the direction of propagation of the original quark. This first 
stage of the process, i.e., the parton cascade development, described by perturbative QCD, 
is terminated at a cut-off value, (A^)cut-off ~ 1 GeV, of the typical transverse momentum. 
Thereafter, the second stage, involving the non-perturbative "confinement" process, takes 
over binding the quarks and gluons into color neutral hadrons. This second stage is usually 
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described by one of the available phenomenological hadronization models such as the LUND 
string fragmentation model |35| or the cluster fragmentation model [36|. Detailed Monte 
Carlo numerical codes now exist |36|, |37], Q which incorporate the two stages outlined above. 



These codes provide a reasonably good description of a variety of relevant experimental data. 
Clearly, however, this is a numerically intensive approach. 



Local Parton-Hadron Duality 

There is an alternative approach that is essentially analytical and has proved very fruitful in 
terms of its ability to describe the gross features of hadronic jet systems, such as the inclusive 
spectra of particles, the particle multiplicities and their correlations, etc., reasonably well. 
This approach is based on the concept of "Local Parton Hadron Duality" (LPHD) ||39|| . 
Basically, in this approach, the second stage involving the non-perturbative hadronization 
process mentioned above is ignored, and the hadron spectrum is taken to be the same, up 
to an overall normalization constant, as the spectrum of partons (i.e., quarks and gluons) 
in the parton cascade after evolving the latter all the way down to a cut-off transverse 
momentum (^l)cut-ofr ~ ~ ^ ew hundred MeV, where R is a typical hadronic size. At 
present the only justification for such an approach seems to be that it gives a remarkably 
good description of the experimental data including recent experimental results from LEP, 
HERA and TEVATRON |[40|| . Justification of LPHD at a more fundamental theoretical 



level, however, is not yet available. Nevertheless, it serves as a good phenomenological tool. 

The main assumption in LPHD is that the actual hadronization process, i.e., the con- 
version of the quarks and gluons in the parton cascade into color neutral hadrons, occurs 
at a low virtuality scale of order of a typical hadron mass independently of the energy of 
the cascade initiating primary quark, and involves only low momentum transfers and local 
color 're-arrangement' which somehow does not drastically alter the form of the momentum 
spectrum of the particles in the parton cascade already determined by the 'hard' (i.e., large 
momentum transfer) perturbative QCD processes. Thus, the non-perturbative hadroniza- 
tion effects are lumped together in an 'unimportant' overall normalization constant which 
can be determined phenomenologically. 

A good quantitative description of the perturbative QCD stage of the parton cascade evo- 
lution is provided by the so-called Modified Leading Logarithmic Approximation (MLLA) 
of QCD, which allows the parton energy spectrum to be expressed analytically in terms of 
functions depending on two free parameters, namely, the effective QCD scale A c g (which 
fixes the effective running QCD coupling strength a° s (Q 2 )) and the transverse momentum 
cut-off Q . For the case Q = A eff , the analytical result simplifies considerably, and one gets 
what is referred to as the "limiting spectrum" [[39, 4^]. For asymptotically high energies 
of interest, i.e., for energies £j et of the original jet-initiating quark satisfying E- ]et 2> A e ff, 
the limiting spectrum can be approximated by a Gaussian in the variable £ = ln(l/x), with 
x = -Eparton/ Ejet, -Eparton being the energy of a quark (parton) in the jet: 
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where Y = In (£ jet / A eff ) , £ « W 2a 2 = (f^) 1/2 with 6 = (11JV C - 2iV»/3 , N c being 
the number of colors and Nf the number of flavors of quarks involved, and N piLrton (Y) ~ 
exp{(16A^ c F/6) 1 / 2 } is the average total multiplicity of the partons in the jet. 

Eq. (HI) gives us the spectrum of the partons in the jet. By LPHD hypothesis, the shape 
of the hadron spectrum, dNh/dx (with x = Eh/E- 3et , Eh being the energy of a hadron in the 
jet), is given by the same form as in Eq. (0), except for an overall normalization constant 
that takes account of the effect of conversion of partons into hadrons. Phenomenologically, 
for given values of A c g and E^ et , the normalization constant can be determined simply from 
overall energy conservation, i.e., from the condition Jq x dNh ^' x ^ dx = 1 . The value of A c g is 
not known a priori, but a fit to the inclusive charged particle spectrum in e + e~ collisions at 
center-of-mass energy E cm = 2£j et ~ 90 GeV (Z-resonance) gives ~ 250 MeV. 

Note that, within the LPHD picture, there is no way of distinguishing between various 
different species of hadrons. Phenomenologically, the experimental data can be fitted by 
using different values of A e fj for different species of particles depending on their masses. For 
our consideration of particles at EHECR energies, all particles under consideration will be 
extremely relativistic, and since, in our case, E }et ~ m x /2 ^> A cff , the hadron spectrum will 
be relatively insensitive to the exact value of A e g. Also, one can safely assume that at the 
asymptotically high energies of our interest, all hadrons — mesons as well as baryons — will 
have the same spectrum. However, the dominant species of particles in terms of their overall 
number will be the light mesons (pions); baryons typically constitute a fraction of < (3 - 
10)% as indicated by existing collider data. For more details on various phenomenological 
aspects of the LPHD hypothesis, see the reviews HO . 



Nucleon, Photon and Neutrino Injection Spectra 



Using the MLLA + LPHD hadron spectrum discussed above, and assuming that each X 
particle on average undergoes A^-body decay (typically AT < 3) to N q quarks (including 
antiquarks) and A"^ leptons (neutrinos and/or charged leptons), so that N = N q + Ng, and 
assuming that the energy mx is shared roughly equally by the A" primary decay products 
of the X, the nucleon injection spectrum, ^-^(Ei,ti), from the decay of all X particles from 
TDs at any time ti can be written as 



dn x (tj 
dh 



■NJ, 



N 



1 



/27T 



cxp 



-In 2 



x 



(3) 



where E^ denotes the energy at injection, is the number of X particle released per unit 
volume per unit time at time ti, /n is the nucleon fraction in the hadronic jet produced by 
a single quark, x = NEi/mx, x* = (NA eS /mxY^ 2 , and A" norm is the normalization constant 
defined by 



Anorm^x) 
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(4) 



An important point about the nucleon injection spectrum given by Eq. (§) is that, un- 
like the spectrum predicted in the standard diffusive shock acceleration theory (see, e.g., 
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Refs. pO| , [378] , [19[]), the injection spectrum in the TD scenario (or, for that matter, in any 
non-acceleration scenario in which the energetic particles arise from decay of massive ele- 
mentary particles), is not, in general, a power-law in energy. Although, in the energy regions 
of our interest, the spectrum (|3]) can be approximated [|7| by power-law segments (oc E^ a ), 
the power-law index a, in the energy regions of interest, is generally smaller than that in 
shock acceleration theories — the latter typically predict a > 2. In other words, the in- 
jection spectrum in non-acceleration theories is generally harder (or flatter) compared to 
that in conventional acceleration theories. This fact has important consequences; it leads 
to the prediction of a pronounced "recovery" |§ of the evolved nucleon spectrum after a 
partial GZK "cut-off" and the consequent flattening of the spectrum above ~ 10 11 GeV. A 
relatively hard spectrum may also naturally give rise to a "gap" in the measured EHECR 
spectrum fH}| . 

The photon injection spectrum from the decay of the neutral pions (tt° — > 2j) in the jets 
is given by 



$ 1 (E i ,t i )~2 —^o(E) , (5) 

JEi hi 

where $„-<) ~ 1 1 ~^ JV $jv is the neutral pion spectrum in the jet. 

Similarly, the neutrino (f At + ^) injection spectrum resulting from the charged pion decay 
[7r ± — > /i ± z/ At (z/ At )] can be written as HH [| 

r mx/N dE 

^^-2.34/ — -$ (7r++w -)(S , (6) 

J2.UEi E 



where $ (7r++7r -) ~ f ^^iv- 

The decay of each muon (from the decay of a charged pion) produces two more neutrinos 
and an electron (or positron): ^ — > e ± v e {p^)v^{y^). Thus each charged pion eventually 
gives rise to three neutrinos: one z/ M , one and one v e (or z/ e ), all of roughly the same 
energy. So the total + injection spectrum will be roughly twice the spectrum given in 
Eq. (^|), while the total v e + v e spectrum will be roughly same as that in Eq. (|6]). 

Note that, if the hadron spectrum in the jet is generally approximated by a power-law in 
energy, then nucleon, photon and neutrino injection spectra will also have the same power- 
law form all with the same power-law index. 

It is worth emphasizing here that while using the LPHD hadron spectrum in the analysis 
of the TD scenario of EHECR one should keep in mind that there is a great deal of uncertainty 
involved in extrapolating the QCD based hadron spectra — which have been tested so far 
only at relatively 'low' energies of ~ 100 GeV — to the extremely high energies of ~ 10 15 GeV 
or so. 

Evolution of the Particle Spectra 
Nucleons 

The evolution of the nucleon spectrum is mainly governed by interactions of the nucleons 
with the CMBR. The relevant interactions are pair production by protons (pjb pe + e~), 
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photoproduction of single or multiple pions by nucleons iV (N^i, —>■ Nmr, n > 1), and 
neutron decay. Here 7^ stands for a background photon, in this case, a CMBR photon. At 
EHECR energies, the photopion production is the dominant process. This is a drastic energy 
loss process for nucleons, and is the basis of the well-known prediction of the GZK cut-off |34 



of the evolved EHECR nucleon spectrum at energies above ~ 10 11 GeV. This process also 
limits the distance of a possible source of the observed EHECR particles to distances less 
than ~ 50Mpc M, 



7-rays 

The 7-rays at EHECR energies interact via pair production (PP: 77;, — > e + e~) and double 
pair production (DPP: 77^ — > e + e _ e + e~), while the electrons (positrons) interact via inverse 
Compton scattering (ICS: e^b — > e'j) and triplet pair production (TPP: e^b — * ee + e~). 
In addition, the electrons (positrons) suffer synchrotron energy loss in the extragalactic 
magnetic field (EGMF). The background photons (7^) involved in the PP process are mainly 
the universal radio background (URB) photons for 7-rays above ~ 10 19 eV, the CMBR 
photons for 7-rays between ~ 10 14 eV and ~ 10 19 eV, and the infrared and optical (IR/O) 
background photons for 7-rays below ~ 10 14 eV. 

The evolution of the 7-ray spectrum is complicated due to the fact that PP and ICS 
processes together lead to development of electromagnetic (EM) cascades, whereby any elec- 
tromagnetic energy in the form of 7-rays and/or electrons (positrons) released at an energy 
which is above the threshold for PP process on photons of a particular background, cascades 
down to progressively lower energies through a cycle of PP interactions (on the photons of 
progressively higher energy backgrounds) and ICS interactions (mainly on the CMBR). The 
cascading has the overall effect of increasing the 7-ray flux at EHECR energies because it 
causes an effective increase of the attenuation length of these 7-rays. Further cascading by 
any cascade photon stops either if the remaining path length of the cascade photon, as it 
propagates from its point of creation to the observation point, becomes less than the mean 
free path for the PP process on the relevant background photons, or if the energy of the 
propagating cascade photon falls below the threshold for the PP process. Thus, depend- 
ing on the distance at which they are first injected, 7-rays of EHECR energies can, after 
propagation, give rise to a 7-ray spectrum that may span the energy range from (say) a few 
tens of MeV all the way up to EHECR energies. This, of course, means that any model 
of electromagnetic energy injection at EHECR energies has to meet the constraint that the 
resulting cascade 7-ray flux at lower energies should not exceed the measured flux at those 
energies. In the context of TD models of EHECR, this was first pointed out in Ref . |2^| . 

The mean attenuation length of EHE 7-rays depends strongly on the density of URB 
and on the strength of the EGMF, both of which are uncertain at the present time. The 
EGMF typically inhibits cascade development because of synchrotron cooling of the e~e + 
pairs produced in the PP process. Depending on the strength of the EGMF, the synchrotron 
cooling time scale may be shorter than the time scale of ICS, in which case the e~ or the e + 
under consideration loses energy through synchrotron radiation before it can undergo ICS, 
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and thus cascade development stops. In this case the 7-ray flux is determined mainly by 
the "direct" 7-rays, i.e., the ones that originate at distances less than the absorption length 
due to PP process. The energy lost by synchrotron cooling does not, however, disappear - 
rather, it appears at a lower energy and can initiate fresh EM cascades by interacting with 
the photons of a higher energy background such as CMBR or IR/O depending on its energy. 
So the overall effect of a relatively strong EGMF is to deplete the 7-ray flux above some 
energy in the EHE region and increase the 7-ray flux below a corresponding energy in the 
'low' (MeV - GeV) energy region, where it will have to meet the constraints imposed by the 
measured extragalactic diffuse 7-ray background ^3|, |44 |. 



In addition to uncertainties in the strength of EGMF and the URB, another major source 
of uncertainty in the cascade calculation is the poorly known IR/O background (see, e.g., 



Ref. [[f5[] for a recent discussion of IR/O backgrounds). The latter strongly influences the 
cascade spectrum in the energy range from ~ 10 11 eV to ~ 10 14 eV. Below this range, 
however, the cascade spectrum becomes relatively insensitive to the model parameters 



that determine the IR/O background. This is a fortunate circumstance because this implies 
that the constraints on TD models derived [jn| [H], [I7|, [5| by comparing the measured 100 
MeV - 10 GeV 7-ray background with TD model predictions are relatively insensitive to 
uncertainties in our precise knowledge of the IR/O background, and hence are fairly robust. 

Neutrinos 

EHE neutrinos suffer absorption through fermion-antifermion pair production on the thermal 
background neutrinos [y + — > f f), where / = e, /x, r, u, q, and Vb is a thermal background 
neutrino. Due to this absorption process, neutrinos of present observed energy E u would 
have to have been injected at redshifts less than z a (E U:0 ) ~ 3.5 x 10 2 (10 20 eV/i?^o) 2 ^ 7 , for 
£„, >3x 10 14 eV !. 

Actually, the /i's, r's, and quarks created in the above absorption process generate further 
neutrinos through decay of the /i's and r's and through decay of the charged pions created 
by quark fragmentation. This leads to a "neutrino cascade" p7j, effectively increasing the 



size of the "neutrino horizon" of the Universe. This, in turn, has the effect of increasing 
the overall neutrino flux around 10 20 eV by a factor of few relative to the case when the 
cascading effect is not taken into account. 

Recently, it has been pointed out j|7| that if neutrinos have a small mass m„ in the eV 
range, then EHE neutrinos (antineutrinos) of energy ~ 4 ( eV/m u ) x 10 12 GeV will annihilate 
on the relic antineutrinos (neutrinos) to produce the Z-boson with a resonant cross section 
of ~ 10~ 32 cm 2 . The hadronic decay of the Z will produce additional photons, neutrinos and 
nucleons, which will add to the photon and neutrino cascading processes mentioned above. 
This process may also leave specific signatures on the EHE neutrino spectrum, the detection 
of which may in the end provide an indirect signature of the neutrino mass and hence of dark 
matter. Detailed self-consistent calculation of nucleon, photon and neutrino fluxes, however, 
remain to be done. 
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Predicted Particle Fluxes and Constraints on TD models 



As discussed above, the predicted particle fluxes depend on a number of parameters: the 
X particle mass nix, the strength of EGMF, the URB, the injection spectra, and so on. 
Recently, detailed numerical calculation of the particle spectra, especially the spectrum of 7- 
rays in the entire energy range from ~ 10 s eV to ~ 10 25 eV, have been done [|15|, 16, 17| taking 
into account the effects of electromagnetic cascading and EGMF. Fig. 1 shows the predicted 




Figure 1: Predicted fluxes of 7-rays (dash - triple dotted line), protons and neutrons (dash- 
dotted lines) and {y^ + u^) (solid lines) for TD models with p = 1, m x = 10 16 GeV and 
EGMF of 10 -12 Gauss. The neutrino curve marked "SLSC12" corresponds to the calculation 
of Ref. while the curves marked "BHSO" and "BHS1" represent the neutrino fluxes from 
Ref. and correspond respectively to p = and the p = 1 models. Also shown are the 
estimates of the atmospheric neutrino background for different zenith angles (hatched 
region marked "atmospheric"). Data points with error bars represent the combined cosmic 
ray data from the Fly's Eye and AGASA experiments [T2| above 10 19 eV, and the thick 



solid line represents piecewise power-law fit to the observed charged CR flux. The dash 
- triple dotted line on the left margin represents experimental upper limits on the diffuse 
7-ray flux at 100 MeV - 5 GeV from EGRET data [^, Points with arrows pointing 
downwards represent approximate upper limits on the diffuse neutrino flux from Frejus [[49 



the EAS-TOP ||50|| , and the Fly's Eye |[T| experiments, as indicated. The projected limit 
shown for the proposed Auger experiment assumes the acceptance estimated in Ref. |52|] for 
non-detection over a five year period. (Courtesy G. Sigl) 
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diffuse particle fluxes for a representative set of values of various parameters involved. Here 
I discuss only the predicted diffuse fluxes of particles assuming uniform distribution of TD 
sources in the Universe. (Possible individual "bursting" sources of EHECR and the possible 
reconstruction of their "images" from data on arrival direction, arrival time and energy 
associated with individual events, leading to information about the source characteristics 
and, in particular about the EGMF, are discussed in the talk by G. Sigl; see this volume.) 

Because the magnitude of X particle production rate is not known a priori, the best one 
can do at the present time is adopt a suitable normalization procedure for the absolute flux 
so as to be able to explain the EHECR data, and then check whether this normalization is 
consistent with all relevant observational data, not just on EHECR particles, but also on 
other relevant particle fluxes at lower energies, especially the diffuse gamma ray background 
measurements in the MeV - GeV region. The normalization of the absolute fluxes in Fig. 1 
is optimized in the sense that it has been determined by fitting the 'observable' particle (i.e., 
the combined nucleon and 7-ray) fluxes to the measured EHECR data by the maximum 
likelihood method |14|] and corresponds to a likelihood significance of > 50%. 

It is clear that TD models, while potentially contributing dominantly to the particle 
fluxes above ~ 10 20 eV, make negligible contribution to cosmic ray flux below ~ 10 19 eV 
because of the relatively hard nature of the particle spectra in these models. Thus the flux 
below ~ 2 x 10 19 eV is presumably due to a conventional acceleration scenario, and was not 
included in the fitting procedure. 

Since pions are the most numerous particles in the jets, their decay products, i.e., photons 
and neutrinos dominate the number of particles at production. Since neutrinos suffer little 
attenuation and can come to us unattenuated from large cosmological distances (except for 
absorption due to e + e~ pair production by interacting with the cosmic thermal neutrino 
background, the path length for which is ^> 100 Mpc), their fluxes, as expected, are the 
largest among all particles at the highest energies. However, their detection probability 
is much lower compared to those for protons and photonsQ. Photons also far outnumber 
nucleons at production, and depending on the level of the URB and EGMF, may dominate 
over the nucleon flux and thus dominate the 'observable' particle flux at EHECR energies. 

An important point to note is that photons and neutrinos in the TD scenario are primary 
particles in that they are produced directly from the decay of the pions in the hadronic jets. In 
contrast, photons and neutrinos in conventional acceleration scenarios can be produced only 
through secondary processes — they are mainly produced by the decay of photoproduced 
pions resulting from the GZK interactions of primary HECR nucleons (produced by the 



acceleration process) with CMBR photons [42|. Of course, these secondary neutrinos and 



photons would also be there in the TD scenario, but their fluxes are sub-dominant to the 
primary ones at the highest energies. 

The shapes of the EHE nucleon- and 7-ray spectra In the TD scenario are "universal" 
|3|, in the sense that they are independent of specific TD process even though different 

4 The EHE neutrinos of TD origin would, however, be potentially detectable by the proposed space-based 
detectors like OWL and AIRWATCH, and ground based detectors like Auger, Telescope Array, and so on; 
see articles on these detector projects in this volume. 



14 



TD sources evolve differently (as parametrized by the parameter p in the X particle pro- 
duction rate). This is because, at these energies, the attenuation lengths of nucleons and 
7-rays are small (< 100 Mpc) compared to Hubble length and so the effects of unknown cos- 
mological evolution of the TD sources are negligible compared to propagation effects. The 
universal shapes of the EHECR nucleon and 7-ray spectra reflect their injection spectra, 
which, as discussed earlier, are determined by QCD. Thus, large statistics measurement of 
the EHECR spectrum by future detectors may give us a probe of new physics, and in partic- 
ular QCD, at energies not currently accessible in laboratory accelerators, provided a TD-like 
non- acceleration scenario of origin of EHECR is correct. 

In contrast to the EHECR nucleons and 7-rays, the predicted 7-ray flux below ~ 10 14 eV 
(the threshold for pair production on the photons of CMBR) and the predicted EHE neutrino 
flux depend on the total energy released integrated over redshift and hence are dependent 
on the specific TD model (i.e., specific value of p). In particular, the 7-ray flux below 
~ 10 11 eV scales as the total electromagnetic energy released from X particles integrated over 



all redshifts and increases with decreasing value of p. This has been used to constrain |53 
TD models from considerations of CMBR distortions and from independent considerations 
of modified light element abundances due to 4 He photodisintegration; for example, this rules 



out |53| the p = TD model. 

The 7-ray flux in Fig. 1 is consistent with the estimates of the upper limit on the back- 
ground in the 100 MeV to ~ 5 GeV from EGRET data f43|. Very recently, new estimates 



of the background flux have been presented which now extend up to ~ 100 GeV with 



roughly the same power-law index as the earlier estimates | 43fl . If the EGMF is signifi- 
cantly larger than the value assumed (~ 10~ 12 Gauss) for the calculation of Fig. 1, then the 
electromagnetic cascade energy transferred to the relevant low energy region will be higher, 
and in this case the 7-ray flux in Fig. 1 will only be marginally consistent, or may even 
be inconsistent, with the experimental diffuse background flux estimates. Note, however, 
that the fluxes in Fig. 1 correspond to the case mj = 10 16 GeV. Lower values of mx are 
possible, which will give lower contribution to the low energy diffuse flux while at the same 
time producing enough energy in the form of X particles to explain the observed EHECR 
flux. This can be seen as follows: 

As we shall discuss in more details in the next section, for a relatively hard power-law 
photon injection spectrum cx E~ a with a < 2, the energy injection rate Qq required to match 
a given differential EHECR flux at a given energy decreases with decreasing value of mx 
(provided, of course, m x > 10 11 GeV, the energy of the highest energy event). For example, 

it is easy to see (see next section) that for a photon injection spectrum oc E~ 1,5 , the value 

1/2 

of Qq required to explain the EHECR flux is roughly proportional to m^ ■ And since the 
cascade 7-ray flux in the < 100 GeV region is essentially directly proportional to Qo injected 
at EHECR energies, a reduction by a factor of 10, for example, of the predicted 7-ray flux in 
the < 100 GeV region is easily achieved by reducing mx by about two orders of magnitude, 
i.e., for mx 10 14 GeV. And, of course, as already mentioned, the low energy 7-ray flux is 
also reduced if the EGMF strength is lower. 

Thus, it seems that TD models of EHECR with m x ~ 10 16 GeV and/or high EGMF 
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(~ 10~ 9 Gauss) are somewhat difficult to simultaneously reconcile with EHECR data and 
the low energy diffuse 7-ray background. On the other hand, models with mx ^ 10 14 GeV 
and low EGMF (< 1CT 11 Gauss) can explain the EHECR data and at the same time be 
consistent with all existing data. 

The neutrino flux indicated by curve marked "SLSC12" in Figure 1 corresponds to over- 
all X particle production rate obtained by the maximum likelihood normalization of the 
'observable' (i.e., photon plus nucleon) flux to the EHECR data for the p = 1 TD model. 
The predicted level of the neutrino flux, therefore, depends on the value of EGMF, the ra- 
dio background, etc., which go into determining the photon flux. The earlier flux estimate 
(the curve marked "BHS1") in Fig.l is higher simply because the overall X particle produc- 
tion rate was normalized to a higher value; that normalization was obtained by normalizing 
the predicted proton (as opposed to photon) flux with the measured flux at a lower energy 
~ 5 x 10 19 eV (where the measured cosmic ray flux is higher than that at EHECR energies) 
- the Fly's Eye and AGASA highest energy events beyond 10 20 eV |T2|] were not yet discov- 
ered at that time! The p = model is ruled out by the upper limits from Frejus as well as 
EAS-TOP experiments. (Actually, the p = model is also ruled out from considerations ||5B 



of 4 He photo-disintegration and CMBR distortions as mentioned above — it corresponds to 
unacceptably high rate of energy injection in the early cosmological epochs.) On the other 
hand, the new neutrino flux estimates for the p = 1 model are consistent with all existing 
experimental upper limits. At energies > 10 20 eV, the predicted neutrino flux in the p — 1 
TD model also dominates over the predicted flux from blazars/AGNs as well as over the 
predicted flux of "cosmogenic" neutrinos produced by interactions of UHE cosmic rays with 
CMBR tfpM (not shown in Fig. 1). For more details on the detectability of neutrino fluxes 



in TD models, see, e.g., Ref. pL E7J, 55, 56 



X-Particle Production Rate Required to Explain the 
EHECR Flux: A Benchmark Calculation 

To have an idea of the kind of numbers for the required X particle production rate, we can 
perform a simple (albeit crude) benchmark calculation as follows: 

Since in TD models, photons are expected to dominate the observable EHECR flux, let 
us assume for simplicity that the highest energy events are due to photons. Let us assume 
a typical three-body decay mode of the X into a qq pair and a lepton: X — > qq£. The two 
quarks will produce two hadronic jets. Let /„- denote the total pion fraction in a jet. Then 
the photons from the two jets carry a total energy -E 7i Totai — (| x 0-9 x §) m x(/7r/0.9) = 
0.2mx(/7r/0.9). Let us assume a power-law hadronic fragmentation spectrum with index 1.5. 
Then the photon injection spectrum from a single X particle can be written as 



diV 7 3 f3E, 



-1.5 



x0.3R P , (7) 
dE y m x \m x J \0.9J ' v ; 

which is properly normalized with the total photon energy -E-vTotai- We shall neglect cosmo- 
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logical evolution effects since photons of EHECR energies have a cosmologically negligible 
path length of only ~ few tens of Mpc for absorption through pair production on the radio 
background. 

With these assumptions, the photon flux j 7 (-E' 7 ) at the observed energy E 1 is simply 
given by 

i i 1 , . „ „ dn x dN~ , , 

h(El ^-ME,)^—^ (8) 

where A(_E 7 ) is the pair production absorption path length of a photon of energy E^. 

Noting that dnx/dt = Qo/mx, and normalizing the above flux to the measured EHECR 
flux corresponding to the highest energy event at ~ 3 x 10 20 eV, given by j(3 x 10 20 eV) ~ 
5.6 x 10~ 41 cm -2 eV -1 sec -1 sr _1 , we get 



Qcrequixed « 2.1 x 10" 21 eV cm" 3 sec^ 1 — — ^ — ^ - , (9) 



10Mpc\ / m x ^ 1/2 
N A 7i3 oo J UO 16 GeV, 

or 

„ 2.1 x 10- cm-sec- (^) ( SI Y 1 ' 2 , (io) 

where A 7i 3oo is the absorption path length of a 300 EeV photon (1 EeV = 10 18 eV). The 
subscript stands for the present epoch. 

The above numbers are probably uncertain by up to an order of magnitude depending 
on the exact form of the injection spectrum, the absorption path length of EHECR photons, 
electromagnetic cascading effect, and so on. Indeed, since electromagnetic cascading effect 
(which we have neglected here) generally increases the photon flux, the above numbers are 
most likely overestimates. Nevertheless, they do serve as crude benchmark numbers. These 
numbers indicate that in order for TDs to explain the EHECR events, the X particles must 
be produced in the present epoch at a rate of ~ 2 x 10 35 Mpc -3 yr _1 , or in more "down- 
to-earth" units, about ~ 23Au~ 3 yr _1 , i.e., about 20 X particles within a solar system-size 
volume per year. 

In the next section I discuss three plausible specific TD processes and examine their 
efficacies with regard to X particle production and EHECR keeping in mind the above rough 
estimates of the required X particle production rate. 

X Particle Production from Cosmic Strings, Monopoles, 
and Necklaces 



Cosmic Strings 

Let us first recall the salient features of evolution of cosmic strings in the Universe; for a 
review, see Refs. 0. Immediately after their formation, the strings would be in a random 
tangled configuration. One can define a characteristic length scale, £ s , of the initial string 
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configuration in terms of the overall mass-energy density, p s , of strings simply through the 
relation 

p s = v/e s , (ii) 

where /i denotes the string mass (energy) per unit length. 

Initially, the strings find themselves in a dense medium, so they move under a strong 
frictional damping force. The damping remains significant until the temperature falls to 
T < (G/i) 1 / 2 ?], where G is Newton's constant and r\ is the symmetry-breaking scale at which 
strings were formed. [Recall, for GUT scale cosmic strings, for example, rj ~ 10 16 GeV, 
H ~ rj 2 ~ (10 16 GeV) 2 , and Gfi ~ 1CU 6 .] In the friction dominated epoch, a curved string 
segment of radius of curvature r quickly achieves a terminal velocity oc 1/r. The small 
scale irregularities on the strings are, therefore, quickly smoothed out. As a result, the 
strings are straightened out and their total length shortened. The energy density in strings, 
therefore, decreases with time. This means that the characteristic length scale £ s describing 
the string configuration increases with time as the Universe expands. Eventually £ s becomes 
comparable to the causal horizon distance ~ t. At about this time, the ambient density of 
the Universe becomes dilute enough that damping becomes unimportant so that the strings 
start moving relativistically. 

Beyond this point, there are two possibilities. Causality prevents the length scale £ s from 
growing faster than the horizon length. So, either (a) £ s keeps up with the horizon length, 
i.e., £ 8 /t becomes a constant, or (b) £ s increases less rapidly than t. 

Let us consider the second possibility first. In this case, the string density falls less 
rapidly than t~ 2 . On the other hand, we know that the radiation density in the radiation- 
dominated epoch as well as matter density in the matter-dominated epoch both scale as 
t~ 2 . Clearly, therefore, the strings would come to dominate the density of the Universe at 
some point of time. It can be shown that this would happen quite early in the history of 
the Universe unless the strings are very light, much lighter than the GUT scale strings. A 
string dominated early Universe would be unacceptably inhomogeneous conflicting with the 
observed Universe^. 

The other possibility, which goes by the name of "scaling" hypothesis, seems to be more 
probable, as suggested by detailed numerical as well as analytical studies 0, The nu- 
merical simulations generally find that the string density does reach the scaling regime given 
by p s ,scaiin g oc 1/t 2 , and then continues to be in this regime. It is, however, clear that in order 
for this to happen, strings must lose energy in some form at a certain rate. This is because, 
in absence of any energy loss, the string configuration would only be conformally stretched 
by the expansion of the Universe on scales larger than the horizon so that £ s would only 
scale as the scale factor oc t 1 ^ 2 in the radiation dominated Universe, and cx t 2 ^ 3 in the matter 
dominated Universe. In both cases, this would fail to keep the string density in the scaling 
regime, leading back to string domination. In order for the string density to be maintained 



5 However, a string dominated recent Universe — dominated by 'light' strings formed at a phase transition 
at about the electroweak symmetry breaking scale — is possible. Such a string dominated recent Universe 
may even have some desirable cosmological properties ]5^ |. Such light strings are, however, not of interest 
to us in this discussion. 
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in the scaling regime, energy must be lost by the string configuration per unit proper volume 
at a rate p Sj i ss satisfying the equation 

R 

As.total = ~^~ftP s + /Woss , (12) 

where the first term on the right hand side is due to expansion of the Universe, R being the 
scale factor of the expanding Universe. In the scaling regime p s ,totai — — 2p s /£, which gives 
As,ioss = —Ps/t in the radiation dominated Universe, and p Sj i oss = — (2/3)p s /t in the matter 
dominated Universe. 

The important question is, in what form does the string configuration lose its energy 
so as to maintain itself in the scaling regime? One possible mechanism of energy loss from 
strings is formation of closed loops. Occasionally, a segment of string might self-intersect by 
curling up on itself. The intersecting segments may intercommute, i.e., "exchange partners" 
, leading to formation of a closed loop which pinches off the string. The closed loop would 
then oscillate and lose energy by emitting gravitational radiation and eventually disappear. 
It can be shown that this is indeed an efficient mechanism of extracting energy from strings 
and transferring it to other forms. The string energy loss rate estimated above indicates 
that scaling could be maintained by roughly of order one closed loop of horizon size (~ t) 
formed in a horizon size volume (~ t 3 ) in one hubble expansion time (~ t) at any time t. In 
principle, as far as energetics is concerned, one can have the same effect if, instead of one or 
few large loops, a large number of smaller loops are formed. Which one may actually happen 
depends on the detailed dynamics of string evolution, and can only be decided by means of 
numerical simulations. 

Early numerical simulations seemed to support the large (i.e., ~ horizon size) loop for- 
mation picture. Subsequent simulations with improved resolution, however, found a lot of 
small-scale structure on strings, the latter presumably being due to kinks left on the strings 
after each crossing and intercommuting of string segments. Consequently, loops formed were 
found to be much smaller in size than horizon size and correspondingly larger in number. 
Further simulations showed that the loops tended to be formed predominantly on the scale 
of the cut-off length imposed for reasonable resolution of even the smallest size loops allowed 
by the given resolution scale of the simulation. It is, however, generally thought that the 
small-scale structure cannot continue to build up indefinitely, because the back-reaction of 
the kinky string's own gravitational field would eventually stabilize the small-scale structure 
at a scale of order G/jbt. The loops would, therefore, be expected to be formed predomi- 
nantly of size ~ Gfit, at any time t. Although much smaller than the horizon size, these 
loops would still be of 'macroscopic' size, much larger than the microscopic string width 
scale (~ i]^ 1 ~ jU -1 / 2 ). These loops would, therefore, also oscillate and eventually disappear 
by emitting gravitational radiation. Thus, according to above picture, the dominant mecha- 
nism of energy loss from strings responsible for maintaining the string density in the scaling 
regime would be formation of macroscopic-size (^> r/ _1 ) loops and emission of gravitational 
radiation by these loops. More details of this picture of cosmic string evolution can be found 
in Refs. @. 



19 



How is the above picture of cosmic string evolution relevant for cosmic rays? How does 
one get X particles from cosmic strings? 

One way of getting X particles from cosmic strings is through the so-called cusp-evaporation 
mechanism 0. I will not discuss this mechanism here, but the resulting X particle production 
rate generally turns out to be too low to be of relevance to EHECR. 

Another possibility arises as follows: As the closed loops oscillate and emit gravitational 
radiation, they lose energy and shrink. Eventually, when a loop's radius shrinks to a size of 
order the width of the string, the string unwinds and turns into X particles, which will then 
decay, producing high energy particles. However, each loop in the end only produces of the 
order of one X particle. The resulting cosmic ray flux is again too low to be observable ||. 

Clearly, the only way cosmic strings may produce large number of X particles is if macro- 
scopic lengths of string are involved in the X particle production process. One such mecha- 
nism was suggested in Refs. |7], [| based on the following arguments. The picture of gravi- 
tational radiation being the dominant energy loss mechanism for cosmic strings rests on the 
assumption that the loops themselves do not self-intersect frequently. Since the motion of 
a freely oscillating loop is periodic (with a period of L/2, L being the invariant length of 
the loop J58|) , a loop formed in a self- intersecting configuration will undergo self-intersection 
within its first period of oscillation. If the loop does indeed self-intersect and break up into 
two smaller loops, and if the daughter loops again self-intersect breaking up into two even 
smaller loops, and so on, then one can see that a single initially large loop of size L can break 
up into a debris of tiny loops of size 7/ _1 (thereby turning into X particles) within a time scale 
of ~ L. Since the largest loops are expected to be of size < t, one sees that the above time 
scale can be much smaller than one hubble time. In other words, one large loop can break up 
into a large number of tiny loops (X particles) within one hubble time. Such self-intersecting 
loops would not radiate much energy gravitationally because that would require many peri- 
ods of oscillation. In other words, such repeatedly self-intersecting loops would be a channel 
through which energy contained in macroscopic lengths of string could go into X particles 
instead of into gravitational radiation. As shown in Ref. p)8|| , some non-circular loops could 
also be in configurations which would collapse completely into double-line configurations and 
subsequently annihilate into X particles It has also been argued in Ref. pO] that the 



self-intersection probability of a loop increases exponentially with the number of small scale 
kinks on a loop. 

At the time of work of Ref. ||, however, it was not clear as to what fraction of the string 
energy might go into X particles through processes discussed above. Treating this fraction 
to be a free parameter /, i.e., assuming that a fraction / of the energy extracted from the 
long strings per unit volume per unit time went into X particles, the authors of Ref. |§ found 
that the fraction / had to be rather small, / < 7 x 10 -4 ; otherwise, the predicted cosmic ray 
flux from GUT scale cosmic strings would exceed the observed flux. Results of more recent 
calculation of the predicted observable particle flux |T6] (see Fig. 1) correspond to an upper 
limit on / which is about two orders of magnitude lowerf]. 

6 This is due to the fact that the 'observable' particle flux now includes the gamma ray flux in addition 
to the protons — in contrast to only the much lower proton flux considered in Ref. M — and also because 
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The actual value of / is still unknown. But if gravitational radiation, and not massive 
particle radiation, is the dominant energy loss mechanism for cosmic strings, then the fraction 
/ may actually be much smaller than the above upper limit, in which case the flux of cosmic 
rays produced by cosmic strings through X particles would be small and below the observed 
flux. However, rapid conversion of a significant fraction of the energy in macroscopically large 
loops into tiny loops (X particles) through repeated loop self-intersections due, for example, 
to the presence of large number of kinks on the loops |60], and consequent production of 
significant EHECR flux, cannot be ruled out. 

Very recently, the notion gravitational radiation as the dominant energy loss mecha- 
nism for cosmic strings has been questioned by the results of new numerical simulations 
of evolution of cosmic strings [ 33| . Authors of Ref. [|33] claim that if loop production is 



not artificially restricted by imposing a cutoff length for loop size in the simulation, then 
loops tend to be produced dominantly on the smallest allowed length scale in the problem, 
namely, on the scale of the width of the string. Such small loops promptly collapse into X 
particles. In other words, there is essentially no loop production at all — the string energy 
density is maintained in the scaling regime by energy loss from strings predominantly in 
the form of direct X particle emission, rather than by formation of large loops and their 
subsequent gravitational radiation. This new result, while subject to confirmation by inde- 
pendent simulations, obviously has important implications for EHECR. Indeed, in this case, 
the upper limit on the fraction / mentioned above implies severe constraint on GUT-scale 
cosmic strings. From the results of Ref. |[33|| , the string energy density in the scaling regime 
is p SiSC aiing — /V(0.3t) 2 . With X particle production as the dominant energy loss mechanism, 
we immediately see from Eq. ( |T2"D that rate of production of X particles from strings per 
unit volume must be drix/dt ~ 7A(fi/m x )t' 3 in the matter dominated era. Taking, for 
cosmic strings, [i ~ irr] 2 , and taking mx ~ 0.7r/, where 7] is the symmetry-breaking scale, we 
get by using the constraint imposed by Eq. (|T0|), r\ < 4.2 x 10 13 GeV. Thus the GUT scale 



cosmic strings with 77 ~ 10 16 GeV are ruled out by EHECR data if the results of Ref. |33 
are correct. At the same time X particles from cosmic strings formed at a phase transition 
with i] ~ 10 13 - 10 14 GeV are able to explain the EHECR data. Cosmic strings may thus be 
a "natural" source of extremely high energy cosmic rays if massive particle radiation, and 
not gravitational radiation, is indeed their dominant energy loss mechanism. 

Cosmic string formation with rj ~ 10 14 GeV rather than at the GUT scale of ~ 10 16 GeV 
is not hard to envisage. For example, the symmetry breaking scheme SO (10) — > SU(3) x 
SU(2) x U(l) Y x U(l) can take place at the GUT unification scale M GUT ~ 10 16 GeV; with 
no U(l) subgroup broken, this phase transition produces no strings. However, the second 
U(l) can be subsequently broken with a second phase transition at a scale ~ 10 14 GeV to 
yield the cosmic strings relevant for EHECR. Note that these strings would be too light to 
be relevant for structure formation in the Universe and their signature on the CMBR sky 
would also be too weak to be detectable. Instead, the extremely high energy end of the 
cosmic ray spectrum may offer a probing ground for signatures of these cosmic strings. 



of the maximum likelihood normalization of the predicted flux to the observed data. 
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Monopoles 

If monopoles were formed at a phase transition in the early Universe, then, as Hill || 
suggested in 1983, a metastable monopole-antimonopole bound state — "monopolonium" 
- is possible. At any temperature T, monopolonia would be formed with binding energy 
E b > T. The initial radius of a monopolonium would be ~ ^g^/ E b , where g m is the 
magnetic charge (which is related to the electronic charge e through the Dirac quantization 
condition eg m = N/2, N being the monopole's winding number). Classically, of course, 
the monopolonium is unstable. Quantum mechanically, the monopolonium can exist only 
in certain "stationary" states characterized by the principal quantum number n given by 
r = n^a^, where n is a positive integer, r is the instantaneous radius, and = 8a e /rriM 
is the "magnetic" Bohr radius of the monopolonium. Here a e = 1/137 is the "electric" 
fine-structure constant, and m,u is the mass of a monopole. Since the Bohr radius of a 
monopolonium is much less than the Compton wavelength (size) of a monopole, i.e., -C 
rn^j , the monopolonium does not exist in the ground (n = 1) state, because the monopole 
and the antimonopole would be overlapping, and so would annihilate each other. However, 
a monopolonium would initially be formed with n ^> 1; it would then undergo a series of 
transitions through a series of tighter and tighter bound states by emitting initially photons 
and subsequently gluons, Z bosons, and finally the GUT X bosons. Eventually, the cores 
of the monopole and the antimonopole would overlap, at which point the monopolonium 
would annihilate into X particles. Hill showed that the life time of a monopolonium is 
proportional to the cube of its initial radius. Depending on the epoch of formation, some 
of the monopolonia formed in the early Universe could be surviving in the Universe today, 
and some would have collapsed in the recent epochs. It can be shown flO| that monopolonia 



collapsing in the present epoch would have been formed in the early Universe at around the 
epoch of primordial nucleosynthesis. 

The X particles produced by the collapsing monopolonia may give rise to EHECR. This 



possibility was studied in details in Ref. [10], who showed that this process, like cosmic 



strings, can also be described by an equation for X particle production rate described by 
Eq. (|l|) with p = 1. The efficacy of the process, however, depends on two parameters, 
namely, (a) the monopolonium-to-monopole fraction at formation (£y) and (b) the monopole 
abundance. The latter is unknown, while the former is in principle calculable by using the 
classical Saha ionization formalism. However, phenomenologically, since a monopole mass 
can be typically vtlm ~ 40mx (so that each monopolonium collapse can release ~ 80 X 
particles), we see from Eq. (|T0|) that one requires roughly (only!) about 3 monopolonium 
collapse per decade within a volume roughly of the size of the solar system. Whether or not 
this can happen depends, as already mentioned, on £j as well as on the monopole abundance, 
the condition being (tt M hl 00 )h m £f ~ 1.7xl0- 8 (m x /10 16 GeV) 1 / 2 ^ Mpc/A 7i30 o), where 



VIm is the mass density contributed by monopoles in units of closure density of the Universe, 
and hioo is the present hubble constant in units of 100 km sec -1 Mpc^ 1 . Thus, as expected, 
larger the monopole abundance, smaller is the monopolonium fraction £f required to explain 
the EHECR flux. Note that, since must always be less than unity, the above requirements 
can be satisfied as long as (^m^ioo)^ioo > 1-7 x 10 _8 (mx/10 16 GeV) 1 / 2 . Recall, in this 
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context, that the most stringent bound on the monopole abundance is given by the Parkar 
bound [p9| , (fiAf^ioo)parkar ^ 4x 10~ 3 (m M /10 16 GeV) 2 . The estimate of obtained by 
using the Saha ionization formalism [[|, [1(J shows that the resulting requirement on the 
monopole abundance (in order to explain the EHECR flux) is well within the Parkar bound 
mentioned above. The monopolonium collapse, therefore, seems to be an attractive scenario 
in this regard. A detailed numerical simulation of monopolonium formation to determine 
the monopolonium fraction at formation would, however, be useful in this context. 



Necklaces 

A cosmic necklace is a possible hybrid topological defect consisting of a closed cosmic string 
loop with monopole "beads" on them. Such a hybrid defect was first considered by Hind- 
marsh and Kibble |6l ]. Such defects could be formed in a two stage symmetry-breaking 



scheme such as G — > H X U(l) — > H x Z 2 , where Z 2 is the discrete group {-1,1} under 
multiplication. In such a symmetry breaking, monopoles are formed at the first step of the 
symmetry breaking if the group G is semisimple. In the second step, the so-called "Z 2 " 
strings are formed and then each monopole gets attached to two strings, with the monopole 
magnetic flux channeled along the strings. Possible production of massive X particles from 



necklaces has been pointed out in Ref. [11 



The evolution of the necklace system is not well understood. The crucial quantity is the 
dimensionless ratio r = mjvf/(//c?), where tum denotes the monopole mass, \x is the string 
mass per unit length, and d is the average separation between the monopoles. For r< 1, 
the monopoles play a subdominant role, and the evolution of the system is similar to that 
of ordinary cosmic strings. For r > 1, on the other hand, the monopoles determine the 
behavior of the system. The monopoles sitting on the strings tend to make the motion 
of the closed necklaces aperiodic. The authors of Ref. jy] assume that closed necklaces 



undergo frequent self-intersections, leading to monopole - antimonopole annihilation and, 
consequently, release of massive X particles. The X particle production rate for the necklace 
system can also be written in the form of Eq. (p]) with p — 1. The efficacy of the process 
depends on free parameters, r, /x and tum- For appropriate choice of the parameters, the 



required EHECR flux can be obtained. For more details see Ref. [[LlJ] and the article by 
Berezinsky in this volume. 



Topological Defects in Supersymmmetric Theories 

Recently, it has been pointed out [^2] that in a wide class of supersymmetric unified theories, 
the higgs bosons associated with the gauge symmetry breaking can be 'light' — of mass of 
order the soft supersymmetry breaking scale ~ TeV — even though the gauge symmetry 
breaking scale (and hence the mass of the gauge boson) itself may be much larger. The 
topological defects in these theories can, therefore, simultaneously be sources of ~ TeV 
mass-scale higgs bosons as well as supermassive (mass up to ~ 10 16 GeV) gauge bosons. 
The decay of the TeV higgs may give a significant contribution to the observed diffuse 7-ray 
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background above a few GeV, while the supermassive gauge boson decay may explain the 
EHECR. For more details, see Ref. [62]. 



Conclusion 

There is no dearth of specific models of X particle producing processes involving topological 
defects. Almost all the "realistic" processes studied so far can be parametrized in the form of 
Eq. ([l]) with p = 1. The spectra of various kinds of particles produced for all these processes 
would essentially be similar to the ones shown in Fig. 1. Different processes might, however, 
contribute different amounts to the total flux. It is difficult to say which particular process 
may contribute most to the observed EHECR flux. However, in this respect, the cosmic 
string scenario seems to be relatively parameter free, especially if the new results of Ref. |33 



are correct. Experimentally, it will be difficult, if not impossible, to tell which specific TD 
process, if at all, is responsible for the EHE cosmic rays. The best one can hope for is that 
proposed future experiments like OWL, Auger, and so on, may be able to tell us whether 
topological defects (or, for that matter, any non- acceleration mechanism in general) or some 
other completely different scenario is involved in producing the observed EHECR. In any 
event, the prospect of being able to look for signatures of new physics with the proposed 
EHECR experiments is certainly exciting, to say the least. 



Note added: Recently it has been pointed out |63|] that if indeed decays of X particles 
from TD processes are responsible for the EHECR, then the same TD processes occurring 
in the early Universe may also have given rise to the observed baryon asymmetry of the 
Universe through CP- and Baryon number violating decays of the X particles. 
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